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Abstract 

We discuss a mesoscopic mechanism of exchange interaction in ferromagnet- 
normal metal-ferromagnet multilayers. We show that in the case when the 
metal's thickness is larger than the electron mean free path, the relative orien- 
tation of magnetizations in the ferromagnets is perpendicular. The exchange 

energy between ferromagnets decays with the metal thickness as a power law. 
Suggested PACS index category: 05.20-y, 82.20-w 
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Both the experiment and the theory of ferromagnet-normal metal-ferromagnet multilay- 
ers have attracted a lot of attention [ 1_6 l An example of such a structure consisting of 
two ferromagnetic films separated by a nonferromagnetic metallic film is shown in Fig.l. In 
the case when the metal thickness L is much smaller than the electron scattering mean free 
path I the sign of the exchange interaction energy between the ferromagnet's magnetizations 
oscillates as a functions of L with a period of order of the Fermi wave length. As a result the 
magnetic structure of the system oscillates between ferromagnetic and antiferromagnetic ori- 
entations of the ferromagnet's magnetizations t 1-6 !. The explanation of this phenomenon is 
based on the fact that the interlayer exchange energy is due to Ruderman-Kittel interaction 
between electron spins in different ferromagnets. 

In the case of low temperatures and at \r — r'\ » I the exchange energy between two 
localized spins < J{r, r') > averaged over the scattering potential configurations decays 
exponentially with \r — r'\ ' 7 l Here r and r' are coordinates of spins and brakets <> stand 
for averaging over realizations of the scattering potential in the metal and the ferromagnets. 
Recent experiments on ferromagnet-metal-ferromagnet multilayers t 8 ' imply, however, that 
the exchange energy between the ferromagnets does not decay exponentially at L ^> I and 
that the equilibrium relative orientation of the ferromagnet's magnetizations is perpendicular 
independently of L. Phenomenologically, this situation can be described by an effective 
energy per unit area 

E = -J°(m l -m 2 ) + B[(m l -m° 2 ) 2 - 1] (1) 

in the case IB 3> \J°\. Here m\ and m° are, averaged over the film's volumes, unit vectors 
parallel to magnetizations of the ferromagnetic films. Indices 1,2 indicate the first and the 
second ferromagnetic film respectively, J° and B(9) are bilinear and biquadratic coupling 
coefficients. In general, B(9) is a smooth function of the angle 6 between and m°. 

In this paper we discuss a theory of this phenomenon. It has been shown in l 9 " 11 ' that 
the exponential decays of the average < J(r, r') > is connected to the fact that it has a 
random sign at large L. The modulus of the exchange interaction decay with L as a power 
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law. 

We can introduce a local exchange energy J(p) between the ferromagnets as an average 
of J(r, r') over a ferromagnet's surface area of order of L 2 . Here p is coordinate along the 
films. We assume that J(p) is small enough and spatial dependence of the magnetizations 
on the scale of order of L can be neglected. According to Slonczewski t 12 - 13 !, biquadratic term 
proportional to B in Eq. 1 can originate from the existence of spatial fluctuations of the sign 
of exchange interaction J(p) = (J) +5J(p) along the layers. The fluctuations of J(p) cause 
fluctuations of directions of magnetizations. Energy associated with spatial fluctuations in 
the magnetization's directions can be represented as 

E(J(p),mi(p)) = - j d 2 pJ{p){m 1 {p) • m 2 (p)) 

+ad [ d 2 p\ dmi ^ ■ dmi ^ + dm M . d m 2(p) l 
J I dp dp dp dp \ 

where the first term corresponds to the interfilms exchange energy, the second term is as- 
sociated with the gradients of magnetizations inside the films, d is the ferromagnetic film's 
thickness and a is a coefficient characterizing the exchange energy value in the ferromagnets. 

In the case when SJ(p) 3> (J) and J(p) has a random sign, the energy E(J, m.j(p)) has 
a minimum at a sample specific realization ra;(p) = mj + 5rrii(p; [5 J]) with ml-Lm® I 12,14 ' 
and 

B.£jO(») (3) 
G = / <Pp{6J(p)SJ(0)) (4) 

Here B$ is a number of order unity t 14 '. 

Let us consider the case when J(p) has random sign due to mesoscopic fluctuations 
of Ruderman-Kittel oscillations inside the metal I 9-11 !. We assume that the ferromagnetic 



film's thickness d > L s = ^ D/u s and L T = > d,L . The latter inequality allows us to 
neglect temperature dependence of B. Here D is the diffusion constant, which is assumed 
to be the same in the ferromagnetic and nonferromagnetic parts of the sample, and uj s is the 
exchange spin splitting energy in the ferromagnets. We will show that in the case L s < L, 



(5) 



while in the case L s > L, 



(6) 



Here 7 and 71 are smooth functions of 9 of order unity and E c = is the Thouless energy. 
Qualitatively, Eqs.5,6 can be understood as follows: In the case \r — r'\ I the random 
oscillations of J(r,r') exhibit a long range sign correlations l 15 '. In the case L s <C L these 
long range correlations should be cut of at a length of the order of L. As a result, the 
fluctuations of the exchange energy averaged over the area of order of L 2 is of order E c ; 
and they are 5-correlated at a distances larger than L. This leads to Eq.5. We think that 
the estimate presented above can be relevant for the experiment ^. In the opposite limit 
L s ^> L the cut off length is L s . The fluctuations of the exchange energy averaged over the 
area of order L 2 S is of order uj s . This leads to Eq.6, which is independent of L. 

To derive the results presented above we describe the exchange energy splitting in ferro- 
magnet with the help of an effective Hamiltonian 



Here H is the Hamiltonian of free electron gas in a random potential U(r), h(r,9) = 
u s m(r, 6) is the effective magnetic filed which is acting only on electron spins, m(r, 6) 
denotes a unit vector parallel to the magnetic moment in ferromagnets in the case when 
the angle between m\ and m\ is 9, and cr = {<J x ,a y ,a z } is the Pauli matrixes vector. 
We assume the following correlation properties of random potential: < U(r) >= and 
< U(r)U(r') >= 2^7^( r — r ')- Here z/ is density of states at fermi level, r is mean free 
scattering time of electrons. 

To get the correlation function (5J(0)SJ(p)) we consider sample specific fluctuations of 
thermodynamic potential £1(9) of the electrons as a function of 9, 



H = Hv + h(r-9)(T 



(7) 



£1(6) = (Q(9)) + 5Q(9). 



(8) 
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Using the identity = Jd 2 p5J(p) we get 

G(6) = <(^) 2 > (9) 

In the case of noninteracting electrons we can express thermodynamic potential as Q = 

J dfiN(fi) , where N(fx) is the number of electrons at given chemical potential //. Then, the 
o 

correlation function of fluctuations of thermodynamic potential has the form 

(5Q(9 1 )5Q(9 2 )) = J dfi 1 dfi 2 (5N(fi 1 ,9 1 )5N( f i 2 ,9 2 )) (10) 
o 

To calculate it we use the usual diagram technique for averaging over configurations of 
disordered potential ' 16 1. Diagrams for correlation function of number of electrons are shown 
in Fig. 2. As a result we have 

(5N( f x 1 ,9 1 )5N( f i 2 ,9 2 )) = -T^ co n Re f d 3 rd 3 r> ' D^[ry^ n )D v ^ \r;u n ) = 

2 „ r d 2 q . oj n 
= n 6 J J2nf \ (EJK, 9 2 ) + Dq 2 + u n + ifc - ^ 2 )) 2 (U) 

where u) n = 2nnT is the Matsubara frequency, n = 1,2.. and a,j3, 7,/i are spin indices. 
Diffusion propagators D]^(r, r'; tu n ) obey the equation 

[—DA + uj n + - H2)]8at&w + i(h(r; 9 1 )cr 1 ^ a( . - h(r; 9 2 )cr ai 5^) D$(r, r'; u n ) = 

= 8(r-r')8^8 aifi (12) 

The second equality in Eq.ll is the representation in terms of eigenvalues of Eq.12. In the 
case of geometry of the system, shown in Fig.l, the eigenvalues are equal to Dq 2 +E m (9i : 9 2 ). 
Here the spectrum E m (9i, 9 2 ) is determined by equation 

d 2 

[-DS^S^— + i(h(r; 9 l )a 1 ^ aii - h(r; 6> 2 )<T Q ^ 7jU )]* m (2;; /i, f ) = E m ^ m (z; 7, a) (13) 

To calculate Eq.ll we use following equalities 

d 2 q 



d r d 2 q If, 1 dlndet(p) 

d^J j2^4mJ c P + Dq 2 + co n + ifa - fi 2 ) dp " 



L2 

= T~R1 — Mdet(ipo)det(-ip )) (14) 

07TJJ dUJ n 
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Here det(p) = Y\{p—E m ) is the spectral determinant of Eq.13, and p = J §(w n + i(fi\ — /i 2 ). 
In expression Eq.14 the integration contour C runs around zeros of det(p). Let us note that 
although Eqs. 10,11 are formally divergent, their contribution to ((^|p) 2 ) is finite. 

Let us consider the case L s <C L,d when results do not depend on uj s . To define boundary 
conditions for Eq.13 it is convenient to introduce operators 



1 ± ( m(z; ^i)<Ti ) ( m(z; 6» 2 )cr 2 



(15) 



Then the boundary conditions for ^ m (z;^,a) are: -^S + ^ m = at z = ±(d + L/2) and 
S-^ m = at L/2 < \z\ < d+ L/2 . As a result, a solution of the eigenvalue problem Eq.13 
gives the following spectral determinant 

1 + cosh(2dp/L)- 



det(ip) 



sinhpsinh(l + 2d/L)p + 1 - cos(0i - 2 ) 



1 + cosh(2 + 2d/L)p - ^1 + cos(#i + 2 ) 



x 



1 + cosh(2dp/L) 



(16) 



Integrating over the chemical potentials, suming over Matsubara frequencies and considering 
case L T = \Jd/T 3> L; d we get 



1 - 



1 d 



sin 29 



2de\l + $(p) 



1-cos 20 



(17) 



where 



Hp) 



1 + cosh 2dp/L 



(18) 



sinhpsinh(l + 2d/V)p 

The case when L s ~^> L can be studied in the same way giving Eq.6. 
Work of B.S. was supported by Division of Material Sciences, U.S. National Science Foun- 
dation under Contract No.DMR-9625370. 



6 



REFERENCES 
[1] M.Baibich at all, Phys.Rev.Lett.61, 2472, 1988 

[2] S.S.P.Parkin, N.More, K.P.Roche, Phys.Rev.Lett.64, 2304, 1990; ibid 66, 2152, 1991. 

[3] P.Levy, S.Zang, A. Pert, Phys.Rev.Lett.65, 1643, 1990. 

[4] Q.Yang at all, Phys.Rev.Lett.72, 3274, 1994. 

[5] B.Heinrich, J.F.Cochran, Adv. in Physics, 42, 523, 1993. 

[6] Y.Yafet J.Appl.Phys. 61, 4058, 1997. 

[7] Pde Gennes, J.Phys.Rad. 23, 230, 1962. 

[8] P. Fuchs, U. Ramsperger, A. Vaterlaus and M. Landolt, Phys. Rev. B 55, 12546, 1997. 

[9] A. Zyuzin, B. Spivak Pis'ma Zh.Eksp.Teor.Fiz.43, 185, 1986(JETP Lett.43,234). 

[10] Stephen,M.J., E. Abrahams, Solid State Commun.65, 1423, 1988. 

[11] B. Spivak, A. Zyuzin, " Mesoscopic Fluctuations of Current Density in Disordered Con- 
ductors" In "Mesoscopic Phenomena in Solids" Ed. B. Altshuler,U P. Lee, R. Webb, 
Modern Problems in Condensed Matter Sciences vol.30, 1991. 

[12] J.C. Slonczewski, J. Appl. Phys. 73, 5957 ,1993. 

[13] J.C. Slonczewski, Phys.Rev.Lett.67, 3172, 1991. 

[14] E.F. Shender and P.C.W. Holdsworth, Phys.Rev.Lett.76, 2583, 1996. 

[15] A.Zyuzin, B. Spivak, Sov.Phys.Uspehi.41. 195. 1998. 

[16] AAAbrikosov, L.P. Gorkov, I.E. Dzyaloshinski, Methods of Quantum Field Theory in 
Statistical Physics, NY, 1963. 

7 



FIGURES 



L/2+d 
L/2 



-L/2 
-L/2-d 



N 



m 



N 



m 



FIG. 1. A schematic picture of the ferromagnet (F)-normal metal (N)- ferromagnet system 
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FIG. 2. Diagrams for calculation of the correlation function {SN(/j,i, 6i)5N(/j,2, 62)) ■ Solid lines 
correspond to electron Green functions and dashed lines correspond to the correlation function of 
the scattering potential < U(r)U(r') >. 
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